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Chapter 1 
INTROnUCTION 

Advanced composite materials are currently being considered for use 
In a number of applications where the dimensional stability of the 
structure Is of primary Importance. Such applications include several 
space structures: coninur.1 cation antennae, telescopes, solar reflectors, 

satellite power systems, and a space operations center. These materials 
are leading candidates due to their high specific stiffness and low 
coefficient of thermal expansion (CTE) in the fiber direction. Graphite 
fibers actually contract when heated (near room temperature) resulting 
In a small, negative CTE in the fiber direction of a unidirectional 
lamina. In contrast, the CTE perpendicular to the fibers is a rather 
large, positive value coiTiparable to that of aluminum. The CTE of a 
composite laminate may be tailored by selecting the appropriate 
orientation of the individual layers of the laminate. 

in the extremely cold space environment in which these materials 
are to be used, large thermal stresses may develop due to the mismatch 
in the CTE of adjacent plies at different orientations. When the in- 
plane normal stress, becomes sufficiently large, transverse cracks 
result. Such cracks run parallel to the fibers in planes approximately 
perpendicular to the midplane of the ply in which they occur, and form 
at quite even intervals (Fig,. 1). 

Although transverse cracking may not significantly alter the struc- 
tural integrity of a laminate, the CTE may be changed enough to render 
the laminate unacceptable for dimensional stability critical appllca- 
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Fig. 1. Transvrerse Cracking in a [0/90^]^ Laminate 
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tions. Therefore, an understanding of transverse cracking effects on 
the laminate CTE is essential for such applications. 

The pjrpose of this investigation was to study the characteristics 
of thermally-induced transverse cracks, both experimentally and analy- 
tically. Specific enphasis was placed on determining the state of 
stress in the vicinity of a transverse crack, the temperature at which 
transverse cracks form and the subsequent crack density present, the 
occurrence of additional transverse cracking under increased thermal 
load and thermal cycling, and the effects of transverse cracks on the 


laminate CTE 


Chapter 2 
LITERATURE REVIEW 

The p''esence of transverse cracking in con>"osite laminates has been 
a topic of interest for several years. While many studies have focused 
on the parameters affecting the development of tran<:vrrse cracking, only 
recently has its effects on laminate behavior received much attention. 
This chapter presents a brief review of studies involving transverse 
cracking in composite laminates. 

2 . 1 Parameters Affecting Transverse Cracking Under Mechanical Loading 
Although often assumed to be a ply failure, transverse cracking has 
been shown by Bailey and Parvizi [I’J to initiate as localized debondi gs 
at the fiber-matrix interface in glass-epoxy laminates. Und^ f 
increasing mechanical load, the debondi ngs linked up to fnnn tiny rracl s 
a few f'ber diameters in length. These in turn linked up to form trans- 
verse cracks. Bailey et al. [2] reported that for graphite-epoxy and 
glass-epcxy cross-ply laminates (i.e., laminates with fibers alternating 
betv/een 0° and 90° orientations), formation of transverse cracks 
occurred instantaneously in thick transverse layers. In thinner trans- 
verse layers, however, tiny cracks formed, but did not propagate through 
the entire layer until higher strain levels were achieved or did not 
propagate at all prior to final failure. Several authors [3, 4, 5] have 
noted that transverse crack density increases with increasing transverse 
layer thickness. Bailey el al. [2] also reported the occurrence of 
interlaminar delaminations at the n''/90° interface of glass-epoxy cross- 
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p1> laminates. Such eJelaminations were found to occur after the onset 
of transverse cracking huf prior to finai fracture. They did not 
observe dolaminations in graphite-epoxy cross-ply laminates. Crossman 
and Wang [5] studied fttS/90^] graphite-epoxy laminates and found that 
transverse cracking was completely suppressed for transverse layer 
cHcknesses of O.PIO in. (0.026 cm) and below (two plies). Interlaminar 
delaminations were the only form of damage present prior to final fail- 
ure for these laminates. For transverse layers thicker than 0.021 in. 
(0.053 cm), or four plies, traTverse cracking occurred first, followed 
by del aminatien. f-arvizi et al. [6] found transverse cracking com- 
plet'^ly suppressec' in glass-epoxy cross-ply laminates with a transverse 
layer thickness under 0.006 in. (0.015 cm), or one ply. 

Reifsnider and Highsmith [/] observed that transverse crack spacing 
In graphite-epoxy laminates decreased with increasing strain until a 
particular characteristic spacing developed at wt'ich no further cracking 
occurred. Kriz et al. [8] concluded that this minimum crack spacing was 
independent of such factors as load hi‘‘tory and environmental effects, 
Nut dependent only upon the lamina and laminate properties as well as 
the stacking sequence, 

2,2 Thermal ly-1 fidiiced Transverse Cracking 

In addition to mechanically applied loads, transverse cracking may 
result from pure thermal loading. Several authors [4, 9-13] have dis- 
cussed the thermal stresses due to airing at elevated temperatures ami 
factors affecting the magnitude of the residual stresses. Spain [9j 
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discussed several factors inTluenciriij thermally-induced cracking Includ- 
ing resin shrinkage during cure, resin coefficient of thermal expansion, 
curing temperature, and ply orientation. Rowles [14] studied the 
effects of transverse cracks on the thermal response of quasi -isotropic 
graphite-epoxy laminates using nechonically cracked specimens. The 
laminate CIE was found to decrease with increasing crack density. A 
twenty-five percent decrease in CTE was reported at the highest crack 
density 52 cracks/in. (2.05 cracks/mn). 

Hills et a1. [4] studied transverse cracking in graphite-polyimlde 
cross-ply and quasi-isctropic laminates under five types of thermal 
loading and room temperature tensile loading. It was determined that 
transverse crack density and the load to initiate transverse cracks are 
functions of the type of loading and the laminate configuration; cross- 
ply laminates were generally found to exhibit higher ersek densities 
than quasi-isotropic laminates under thermal loading. Casswell [15] 
investigated transverse cracking due to thermal loading In graphite- 
epoxy cross-ply laminates. Transverse cracking was present in all 
laminates after the initial cool-down to -250“F (H6K), the crack den- 
sity increasing as the volume percentage of the transverse layer 
decreased from seventy-five to twenty- five percent. Since tne thermal 
stresses in tne transverse layer (at a given thermal load) increase as 
the volume percentage necreases, lic-wevor, this result may not be attrib- 
uted solely to transverse layer thickness. Under a cyclic type of 
thermal loading in wh’f.b Specimens were cooled in 25'F (14K) increnients 
to -250"F (1161'.), D(jt returned to room temperature for examination after 
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each cooling increment, transverse crack density was found to increase 
with increasing thermal load. 

2.3 Thermal Cycling 

Camahort et a1. [16] thermally cycled unidirectional and 
[04/1^^7^70]^ (a low CTE configuration) graphi te-epoxy specimens twenty- 
five times between -322 and 212®F (76 and 373K) to determine the effects 
on mechanical and thermal properties. Cyclic loading resulted in a six 
to thirteen percent decrease In the transverse tensile modulus, a very 
small increase in the transverse tensile strength and a significant 
decrease in the CTE of the [04/i40/+70]^ laminate. The unidirectional 
laminate was affected very little by the cyclic loading. Fahmy and 
Cunningham [17] thermally cycled selected graphite-epoxy laminates 
between -50 and 302‘‘F (223 and 423K), Mechanical and thermal properties 
were evaluated after 10, 100, ICOO, and 500C cycles. Cyclic loading was 
found to reduce the CTE by as much as cne-third of its original value?. 
Degradation of mechanical properties occurred largely in the first ten 
cycles. Berman [18] thermally cycled composite sandwich structures 
composed of Gf 70/934 graphite-epoxy laminates on aluminum honeycomb core 
between -250 and 150®F (116 and 339K). Transverse cracks were first 
reported after 50 cycles and increased in number until 3000 cycles, 
beyond which only widening of existing cracks occurred. Tensile and 
flexural properties remained virtually unchanged after 15,000 cycles. 

The laminate CTE, however, decreased significantly, mainly In the first 
3,000 cycles. 
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Altnougfi all three studies [15-18] investigated the Influence of 
thermal cycling on the mechanical and thermal laminate properties, none 
correlated these results with the amount of transverse cracking pres- 
ent. 

2 . 4 Transverse Crack Spacing Theories 

A transverse crack spacing theory for cross-ply laminates using a 
modified shear-lag analysis was developed by Garrett and Bailey [3] and 
improved by Parvizi and Bailey [19]. In this theory, a linear elastic 
laminate analysis is used to predict the first transverse crack formed 
ui.-er tensile loading. Uhen the transverse ply cracks, the load it 
can led is transferred to the adjacent longitudinal plies at the plane 
of the crack. The theory assumes that load is transferred l^ack into the 
Lransvers? ply via shear stress at an exponentially decreasing rate 
while moving away from the crack. A maximum stress theory is used to 
predict failure. An envelope of permissii'le crack spacings based on a 
random position of the first crack and a specimen of any given length is 
predicted as a function of applied stress. A comparison of experimental 
results with theoretical predictions for graphite-epoxy cross-ply lami- 
nates by Bailey et a1. [2] showed good correlation. Porvizi and Bailey 
[19] reported an even bettor correlation with glass-epoxy cross-ply 
laminates. Neither study included the residual thermal stresses due to 
cooling when predicting crack spacing. Mills et a1. [4] applied the 
same theory to graphite-polyimide cross-ply laminates and extended the 
theory for use with quasi -i sotropic laminates. Residual stresses v/ere 
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predicted with laminate analysis using stress and temperature dependent 
material properties. Residual stress was added to mechanically applied 
stress when plotting the experimental results. The correlation reported 
for quasi-isotropic laminates was somewhat better than for cross-ply 
laminates. Plunkett [20], using data from Ref. [19], reported that the 
modified shear lag analysis predicts crack density well for thicker 
transverse layers e.g., 0,16 in. (0.4 cm), but begins to over predict 
crack density as the transverse layer thickness decreases. 

Stevens and Lupton [21] proposed a theory based on the assumption 
that crack spacing Is inversely proportional to applied stress. The 
correlation with experimental data for glass-epO)y^ cross-ply laminates 
was quite good, Parvizi and Bailey [19] noted a similar correlation, 
also for glass-epoj^y cross-ply laminates, but only over an intennediate 
range of crack spacing. They also noted that the modified shear-lag 
analysis predicts such an inverse proportionality over an intermediate 
range of cr. ck spacing. 

Wang ar d Crossman [22] and Crossman et a1. [23] extended the energy 
release rate theory of Rybicki et al. [24] to predict transverse crack- 
ing. For a given laminate, a nondimensional available strain energy 
release rate function is determined. This function is reported to be 
quite conplicated, being dependent on tne crack location and geometry, 
stacking sequence, lamina properties and thicknesses, and applied loads 
[23]. Both studies used the finite clement method with a crack closure 
procedure to calculate the available strain energy release rate func- 
tion. This procetlure involves extending a crack by a small, finite 
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amount and computing the work required to close the extension introduced 
[22], This work may bt; thought of as that which is aval lable to extend 
a crack. The strain energy release rate required for crack extension is 
determined experimentally by monitoring a slow, stable crack exten- 
sion. A comparison of the available and required strain energy release 
rates determines whether additional cracking will occur at e given value 
of applied stress. Transverse crack spacing Is predicted by placing two 
cracks at a distance most energetically favorable for a given stress. 
Residual curing stresses are accounted for in the theory by adding their 
effects into the predicted available strain energy release rate. Com- 
parisons with experimental results for graphite-epoxy [+25/90235 
[+25/902 j^ laminates showed good correlations [23], 

2.5 Prediction of Transverse Cracking Effects on the GTE 

Bowles [14] used a classical lamination theory analysis with 
reduced transverse stiffness, E 2 , in the cracked layer to predict the 
effect of transverse cracking on the laminate CTE. A finite element 
analysis was used to determine the relationship between crack density 
and reduction in laminate transverse stiffness. Classical lamination 
theory was then used to determine the reduction in of the cracked 
layer needed to produce the same reduction in laminate transverse stiff- 
ness. A master curve, independent of laminate configuration, for the 
reduction in E 2 as a function of crack density was cons.;ructed and used 
in conjunction with a classical lamination theory analyis to predict 
changes in the laminate CTE due to transverse cracks. Predictions were 
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In good agreement with experimental date for quasi-isotropic graphite- 
epoxy laminates. 




Chapter 3 

ANALYTICAL STUDIES 

In this chapter, the characteristics of transverse cracking are 
studied analytically. The study is begun by attempting to predict the 
teinperature at which transverse cracks initially form in a laminate. 
Next, the state of stress in the vicinity of a transverse crack is 
investigated and the Initial equilibrium crack spacing predicted. The 
formation of additional transverse cracking under increasing thermal 
load is then investigated. Finally, the influence of transverse 
cracking on the laminate CTE is studied. Classical lamination theory, a 
modified shear-lag analysis, and a generalized plane strain finite 
element analysis are used In making these determinations. 

A total of ton laminates were selected for analysis and experimen- 
tal correlation. Six were cross-ply laminate configurations; 

[02/90215, [03/9013, [90/0313, [9O2/O2I5. and [9O3/CI3. Tf.e remaining 
four were quasi-isotropic laminate configurations; [O/+45/9OI3, 
[ 0 / 45 / 90 /- 45 L , [ 90 A 45 /G 1 ,, and [ 90 /- 45 / 0/451 , These laminates were 

selected to study the effects of 90 ® layer thickness and adjacent con- 
straints or. both the formation of transverse cracks and the laminate 
CTE. T!ie material chosen for this investigation was T 300 / 52 C 3 graphite- 
epoxy. 


12 




13 


3.1 General Formulation 

3.1.1 Assumptions 

To render the study of transverse cracking under thermal loading 
and Its effects on the laminate CTE tractable, the following assumptions 
are made: 

1. Individual layers are homogeneous, orthotropic materials. 

2. The laminates are infinite plates in the x and y directions 
(Figs. 2 and 3). 

3. Transverse cracks occur only in the 90° layers. 

4. Transverse cracks occur at evenly spaced intervals and extend 
completely through the laminate. 

As will be seen in the following sections, these assumptions allow for a 
generalized plane strain analysis of transverse cracking using a two 
dimensional, repeating unit cell. 

As shown in Figs. 2 and 3, two perpendicular series of cracks may 
be present in different layers of both the cross-ply and quasi-isotropic 
laminates. To study characteristics of one series of transverse cracks, 
the analysis is limited to a region between two cracks (perpendicular to 
the first series) and at a sufficient distance from each such that any 
effects of these cracks may be neglected. A 3-11 analysis would be 
required to study the interacton between the two series of transverse 


cracks. 
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3.1.2 Temperature Dependence of Material Properties 

In an application \ihere the analysis Is performed over a large 
temperature range, temperature dependent material properties should be 
used. In this analysis the mechanical, thermal, and strength properties 
used are second order least squares curve fits to available data on 
T300/5208 graphite-epoxy. Mechanical and strength property data is 
taken from reference [26], in which all tests were performed on speci- 
mens made from the same prepreg material used in the experimental por- 
tion of this study. Data for the transverse Young's modulus, E 2 » at 
-250°F' (116K) is taken from reference [27]. Thermal property data is 
taken from references [28-31]. Temperature dependence is present in all 
mechanical, strength and thermal properties but to various degrees 
(Appendix A). A constant value for the out-of-plane shear modulus, G 23 . 
end the minor Poisson's ratio, from reference [32] is used in the 

generalized plane strain finite element analysis. Coefficients of the 
least squares curve fits as well as the graphical representation of the 
data and material properties are presented in Appendix A. 

3.1.3 In Situ Transverse Strength of 90° Layers 

The transverse strength of a layer In a nulti directional laminate 
(In situ strength) has been shown to b^; appreciably greater than the 
transverse strength determined from a unidirectional tensile coupon. 
Chamis and Sullivan [33] studied angleplied and quasi - 'sotropic lami- 
nates to determine the in situ strength of the various layers. The 
laminates were mechanically loaded to fracture at which point the stress 
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state of the individual layers v/ere predicted by a classical lamination 
theory analysis. These layer stresses were input into a tensor polyno- 
mial failure criteria and the In situ strengths determined. For AS/350; 
graphite-epoxy, the in situ transverse strength was reported to be three 
times greater than the •. niaxial strength, Flaggs and Kura! (.33] studied 
[O 2 / 9 OJS. C±30/90„]5, and l±60/90^2^, n = 1,2, 3, 4, T300/934 grapMte- 
epo)y^ laminates to determine the in situ transverse strength of the 90* 
layers. Specimens were incrementally loaded and examined for transverse 
cracking after each increment. The laminate load at the onset of 
transverse cracking in the 90® layer was determined and the in situ 
strength of the SO® layers predicted using a classical lamination theory 
analysis. Tfie ratios of the in situ strength of the 90° layers to the 
uniaxial strength reported were 2.48 for the [ 02 / 90215 , 2.21 for the 
[+30/50]g, and 1.S6 for the [+60/90]^ laminate. As the thickness of the 
90® layer increased, the in situ strength in the three types of 
laminates was lower •>nd closer to the strength of a unidirectional 90® 
coupon. 

In this study, the effects of the in situ transverse strength of 
the 90® layers on the predicted temperature of initial transverse crack- 
ing are investigated (section 3.2.2). Determinations of the in Situ 
transverse strengths for the thermally loaded specimens studied experi- 
mentally are presented in section 5.3.1. 
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3.1.4 Failure Criteria for Transverse Cracking 

A maximum normal stress failure criteria is use'J to predict the 
occurrence of transverse cracking. Under tnis criterie, a crack is 
assumed to occur when the transverse normal stress, a^, exceeds the in 
situ transverse strength, of the material. The maximum normal 

stress criteria has been applied to the formation of transverse cracks 
by several other authors [3, 8, 19, 20], Upon formation, transverse 
cracks are assumed to span the entire thickness of the 90° layer, 

3.1.5 Stress Free Temperature 

The temperature at which bonding occurs in a laminate cLnng cool 
down from cure is known as the stress free temperature (SFT). The SFT 
is used as a basis for computing the residual thermal stresses during 
cooling. Although often assumed to be the curing temperature, the SFT 
has been chosen as a lower temperature by some authors. Tsai [35] 
proposed using an unsymmetric [+0]- laminate which warps upon cooling to 
determine the SFT, v/hich is taken to be the temperature at which the 
laminate flattens upon reheating. Pagano and Hahn [35] used this tech- 
nique on T300/5208 graphite-epoxy and determined the SFT to be between 
250 and 300°F (394 ar.vi 422K). A value of 250°F was chosen for their 
analysis. Bowles [37], used the same technique on T300/5208 graphite- 
epoxy specimens made from the same prepreg material used in the experi- 
mental portion of this study. The SFT was determined to be near 
350°F. A value of 350°F is therefore chosen for this investigation. 


3.2 In1t1<Ulon of Transverse Cracks 


3.2.1 Formulation 

A classical lamination theory analysis is used to predict the 
Initial temperature at which transverse cracks form under thermal 
load. A brief summary of the stress-strain-temperature relations for an 
arbitrary laminate are given following the method of Hahn and Pagano 
[llj. A more detailed development of classical lamination theory is 
available in reference [38]. 

The standard laminate configuration used is shown In Fig. 4. 
Assuming the superposition of strains, the total strain in the kth 
layer, e, is given as the sum of the mechanical strain, c™, and thermal 
strain, c 

le(T)}'^ « {c'^CD}'' > le^T)}*^ (3.1) 

where {c^(T}}*^ / {a(t)}^dT (3.2) 

SF 

and T<^p = stress free temperature 

T = temperature of interest 

{a(T)| = temperature dependent CTE 

of the kth layer 


The stress strain relations for the kth layer in the laminate coordinate 
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Fig. 4. Laminate Geometry. 
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lo(T))‘‘ CQ(T)]'^{e"’(T)}'' 


(3.3) 


Substituting (3.1) and (3.2) into (3.3) gives 

{a(T)l^ * CQ(T) j’'({n (T)| - / (ct(T)}*^ch) (3.4) 

^SF 

Equation (3.4) shows that the Inclusion of temperature dependent mater- 
ial properties requires a knowledge of the elastic properties, 

•“* k 

CQ(T)] , at the temperature of interest and the thermal properties, 
(a(T)} , as a function of temperature between the stress free 
temperature and the tetnpcrature of interest. Expressing the total 
strain, {e(T))^‘, in terms of the laminate midplane strains, lc^(T)} , and 
curvatures, {<(!)], gives 


{0(T)}*' = [Q(T)]''{1cq(T)}.-. 2^{^(T)} - I le(T)l‘')ck 

'SF 


Integrating through the 

thickness of 

the lamina 

te yie' 


N(T) + n'^(T)| 


’a(t) 

8(T)‘ 

1 



M(T) w mT(T)| 
1 / 


B(T) 

D(T) 


! 1 
►^'(T) j 


(n‘^(T), M^(T)) n [(HT)]*" /’ {a(T)}‘=dx'l,2)d2 

-h 1' r 


(3.5) 


(3.6) 


where 
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([A], [B], CDl) = f*" [Q(T)]^l,z,z^)d 2 
-h 

Equation (3.6) relates the applied Inplane forces and moments to the 
inidplane strains and curvatures. Under pure thermal loading, 

N(T) “ M(T) = 0 and Eqn. (3.6) is Inverted, giving the laminate response 
under a given thermal load 



L 

'a(T) 

B(T) ' 

-1 

N^(T) 

k(T) 


_B(T) 

0(T) 


m"(t) 


The stresses In a particular layer are determined b> substituting the 
values for Eq(T) andtc(T) into Eqn. (3.5). 

3.2.2 Determination of Crack initiation Temperatures 

Since both the stress and uniaxial strength vary with temperature, 
a trial and error procedure is used to determine the temperature of 
initial cracking. A value for the In situ strength is selected in terms 
of a multiple of the uniaxial strength (i.e., = Cf^(T), where C. is 

a constant). A temperature is selected and the stress and in situ 
strength are calculated. If not equal, another temperature is estimated 
and the procedure rep.ated until the temperature at which 
has been predicted to within 1°F. 
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Since the stresses predicted 1 n symmetric lamlnotes with a classi- 
cal lamination theory analysis are independent of stackinp seque.ice, 
only three analyses need to be performed to represent all ten laminates 
of Interest; one for the [O/QOj]^ and [903/0]^ laminates, another for 
the [Oj/QO]^ and [ 90 / 03 ]j laminates, and a third for the [02/902]^, 
[902/0215 and the four quasi-isotrcpic laminates (since 0 ^ in the 90 “ 
layers of a [O2/9O2I5, [902/0215, and I’uasl-isotropic laminates are 
Identical ). 

The predicted temperatures at which transverse cracks Initially 
form in the three groups of laminates are shown In Fig. 5 for four 
values of the in situ transverse strength. As the percentage of the 90 “ 
layer In the cross-ply laminates decreases, the thermal load needed to 
produce transverse cracking at a given In situ strength decreases. This 
result is expected, since the corresponding increase in the percentage 
of 0“ layer produces a greater constraint on the 90 ® layer and therefore 
a nigher transverse tensile stress, o,, in the 90 “ layer at a given 
thermal load. Increasing the In situ strength results in more than a 
linear increase in the thermal load required to produce cracking, a 
result of the temperature dependent material properties. For the case 
where the in situ strength of the 90 “ layer is taken as the uniaxial 
strength of the material (<^2^^^ = Y^(T)), transverse cracking Is predic- 
ted to be present in all laminates at room temperature. However, as 
will be shown in Chapter 5 , this was not the case. All laminates were 
determined to be free of transverse cracks at room temperature, indicat- 
ing an in situ strength greater than the unidirectional strength, Y^(T). 


TEMPERATURE 
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3 . 3 Stresses In the Vicinity of a Transverse Crac k 

Once the thermal load required to produce transverse crocking Is 
determined using classical lamination theory, an understanding of the 
state of stress in the 90° layer In the vicinity of a transverse crack 
Is necessary to predict the equilibrium crack spacing present and the 
additional thermal load required to produce additional transverse crack- 
ing. To study the stress state, a single transverse crack is assumed to 
form In the 90° layer of the laminate at the initial cracking tempera- 
ture (determined by classical lamination theory) as shown In Fig. 6. At 
the plane of the crack (y=0), the transverse normal stress, io the 
90° layer Is zero. Moving away from the crack In the 90° layer 
(increasing y), stress is transferred back Into the cracked 90° layer 
via shear stresses at the interfaces with adjacent layers. At some 
distance from the crack plane, c »2 fP^ches the in situ transverse 
strength of the layer and another crack is formed. This distance v/ill 
be called the initial equilibrium transvere crack spacing. Transverse 
cracks are assumed to firm periodically act oss the 90° layer at this 
spacing (Fig. 2 and 3). Thus, the task of determining the initial 
equilibrium crack spacing becomes one of determining the transverse 
norma! stress in the 90° layer in the vicinity of an existing transverse 
crack as a function of temperature. 

In this study, a generalized plane strain finite element analysis 
is used to determine th.- state of stress. Results are compared to those 
predicted with a inodifKd shear-lag analysis. The <^''°ss-ply 

laminate is chosen for this investigation as being representative and 
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Fig. 6. InUia'I Equilibrium Crack Spacing Model 
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because It was found to be well behaved experimentally (see Chapter 
5). A brief description of the finite element and modified shear-lag 
analysis methods are presented in the following sections. 

3.3.1 Finite Element Analysis 

Stresses and strains are assumed to be independent of the out-of- 
plane X coordinate (Fig. 6) allowing for a generalized plane strain 
finite element analysis. In such an analysis, a constant value of the x 
direction strain, e , Is assumed and the displacement field is of the 

XX 

form 


u(x,y, 2 ) = U(y,z) h x 

v(x,y,z) = V(y,z) (3.9) 

w(x,y,z) = W{y,z) 

where U, V, and 11 are the unknown nodal displacements. In the case of 
pure thermal loading, is the free laminate thermal ’.train due to 
cool down from the SFT and is determi.ned from a classical larriination 
theory analysis of the thermally loaded laminate in the uncracked 
state. Although may not be constant throughout the laminate in the 
presence of transverse cracks, Bowles [25] found that reductions predic- 
ted in due to transverse cracks had a negligible effect on the y-z 
plane behavior and therefor.e the constant, classical lenination theory 
value could be assumed. 
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The finite element program employed was a modified version of a 
generalized plane strain program for laminated composites. Four node 
isoparametric, quadrilateral elements were used to model the region of 
interest, A complete description of the finite element formulation Is 
given by the program author In reference [39j|. 

The modeled region Is chosen to extend a sufficient distance from 
the crack plane to where stresses and strains are not Influenced by the 
presence of the crack and converge to those predicted by classical 
lamination theory. For the [02^^’^2^s distance was deter- 

mined to be 0.20 in. (0.51 cm). Since the laminates considered in this 
study are symimetric, only one-half of the laminate is modeled. 

The displacement boundary conditions for the modeled region 
(Fig. 7) are: 

U(0,z) = 0 (3.10) 

U(8,z) =0 ~ < z < H 

V(0,z) = 0 

V(B,z) = V* (constrained), T < z < H 
W(y,0) = 0 

The U and V boundary conditions at the line y=0 satisfy the Kirchhoff 
assumption that normals to the laminate midplane remain normal to the 
midplane under unifcrm thermal loading. Hie condition U(B,z) = 0, 

T < z < H was verified by modeling a larger portion of the laminate and 
examining the U displacements in the interior of the grid where this 
conoition Is prescribed. A description cf this procedure is given later 
In this section. 
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The condition V(B,z) = V*, T < z < H accounts for the symmetry 
about the transverse crack plane while allowing for free thermal 
strain. With this condition, the line y=B, T < z < H, is allowed to 
translate in the y direction, but is constrained to remain straight and 
vertical, Uie free thermal displacement, V*, is unknov^n and nyst he 
solved for in the analysis. This constraint boundary condition requires 
a modification to the standard finite element formulation. As illustra- 
ted in Appendix B, this constrained displacement formilation may be 
handled by a systematic modification of the global stiffness matrix and 
global force vector resulting in a system of equations with fewer 
unknowns than the unconstrained set. 

The boundary condition W(y,C} = 0 enforces the laminate syfiimetry 
abcut Its midplane v/hen modsling only the top half of the laminate. 

The modeled region used to study the state of stress in the pres- 
ence of a series of transverse cracks at a uniform crack spacing is 
shown in Fig. 8. Since a line of symmetry exists at t 'e inidplane 

between existing c'acks (y=C), only one-half of the crack spacing is 

modeled. It is noted that the boundary conditions on the modeled region 
are the same as those used in determining the initial equilibrium crack 
spacing (Fig. 7) although the reason for employing some are different. 
The condition V(0,z) 0 previously invoked to satisfy Kirchhoff s 

hypothesis now enforces the line of symmetry at the mi '’plane oetween 

cracks. The conditions U(0,z) = 0 and U(B,z) = 0, T < z < H were 

determined by modeling three crack spacings (Fig. 9) ard examining the U 
displacements at these locations in the Interior of the mesh. It was 


31 




^^IGlNAl 

POOR 


P'^Q'z t'3 
QUALrty 




I — modeled region 



repeating -J 
ceil 
(2B.) 


Fig. 8. Modeled Region for Multiple 
Transverse Cracking. 


32 


ORlGlfiAL PAGS IS 
OF POOR QUALITY 



Fig. 9. Three Crack Spacing Mesh for Confirming 
Boundary Conditions (Shown in Deformed 
State for [O^/^O]^ Laminate). 
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determined that although the U displacements are non-jero in the 
interior of the modeled region for the quasi-isotropic lai.i1nates, they 
converge to zero at y=0 and y»B, T < z < H for all laminates 
investigated in this, study. The remaining boundary conditions are 
invoked for the same reasons as in determining the initial equilibrium 
crack spacing. 

The length of the modeled region, B, may be varied to produce 
different crack spaclngs by changing the y-direction mesh scale factor. 
Some caution riust be exercised, however, to ensure that the element 
aspect ratios do not become so large as to produce an inaccurate 
solution. 

A 272 element, 396 node (918 degrees of freedom) finite element 
mesh (Fig. 10) composed of four node isoparametric quadri lateral ele- 
ments Is used to model the region of interest. Four elements are used 
through the thickness of each ply and the element widths varied along 
the length of the mesn to allow for a greater number of elements near 
the transverse crack. 

Results of the finite clement analyses are presented and compared 
to predictions made with a modified snear-lag analysis following a brief 
description of the latter. 

3.3.2 Modified Shear-Lag Analysis 

The modified sho^r-lag analysis has been used by several authors 
[2, 3, 4, 3, 20] to Pi edict transverse crack density as a function of 
the applied load. A i rlef description of the theory is presented here 
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as used In the analysis of transverse cracking under thermal loading. A 
complete development of the theory is presented in reference [4], 

The theory presented is for a [0^/90^,]^ cross-ply laminate as shown 
in Fig. 11. When transverse crack forms, all thermal stresses at the 
crack plane (y=0) are relieved. IXie to the mismatch in coefficients of 
thermal expansion between the 0“ and 90'* layers, stresses develop in the 
laminate via shear stress between the layers. The shear stress is 
assumed to be proportional to the difference in average elastic dis- 
placements of the 0® and 90° layers 


where 0 is the 90° layer half-thickness (Fig. 11) and H is a constant. 
From a force balance on the 90° layer (Fig. 12', 


Substituting (3.13) into (3.11) gives the differentia! equation for o, 

4 . 


T 


H(v 



(3.11) 


D 



dy 


(3.13) 


Differentiating and using the relations 


dv o 
0 


(3.14) 



and 
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Fig. 11. Modified Shear-Lag Analysis Model. 
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dvgo 02 
dv ^2 

where « Young's modulus of the 0° layer 
E 2 = Young's modulus of thf; 90° layer 
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gives 





Substituting the relations 


and 



EjB + E2D 
B + D 


where ~ Yeung's modulus of the cocflosite 

B = 0° layer thickness 


(3.15) 


(3.16) 


(3.17) 


(3.18.1 


into (3.16) gives a second-order differential equation for cf^ of the 
form 


where 



HE^(B+D) 


(3.19) 

(3.20) 


The constant, H, is determined by using the relation 
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dv ^xy 
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(3.21) 


where ^12 “ shear modulus of 90° layer 


Integrating (3.21) from the 0/90 Interface (v = v^) to the laminate 
midplane (v = VgQ) gives 




(3.22) 


Substituting (3.22) into (3.11) and solving for H yields 


H a - 


'12 


T" 


(3.23) 


Substituting (3.23) into (3.20) gives 


❖ = 


^ 1^2 


(3.24) 


Invoking the boundary conditions 


02i = 0 


•y=0 


I 


' 2 | " ®2 
'large y 


u1 1 


(3.25) 


gives a solution to (3.19) of the form 
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02(y) - exp(-<»^'^^y)) (3.26) 

The shear stress at the interface between adjacent layers is found by 
substituting (3.26) into (3.12) 

Txy(y) = Dtj2“^S^/^exp(.o^^^y) (3.27) 

To determine ©2 in the 90“ layer between two existing cracks at a 
given thermal load, the shear stresses due to the two cracks are 
summed. If the spacing betv;een cracks is S, the shear stress is given 

by 


-txy(y) = fo2'''<.'/^[exp(-<^^^lv) - exp(i^^^(y..S)) (3.28) 

wherr 02 ^ is the 90" layer thermal stress at the given thermal load in 
the abssnse of transverse cracks. Substituting (3.20) into (3.12) and 
integrating between 0 and y gives the desired expression for 

02(y) = <^2^^^ exp(-4>^^^S) - exp(-o^^^y) - exp(<.^^^(y-S)) ] (3.29) 

Since 02 is maximum at the midpoint between existing cracks (y = S/2), 
the next crack will form at this point when 02 reaches 02 ^^^^* There- 
fore, setting 02 = 02 '^^^ at y = S/2 in (3.29) and solving for 02 ^ gives 
an expression for needed to produce additional transverse cracking 
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aj = 02“^ ^[1 + exp(-4^'^S) - 2exp(-^^/^ |)]"^ (3.30) 

The thermal load corresponding to may then be determined using a 
classical lamination theory analysis. 

3.3.3 Analysis of [02/902]j Laminate 

As previously mentioned, the [02/90235 laminate is chosen to 
investigate the state of stress about a transverse crack. For this 
analysis, the temperature of initial transverse cracking is taken as the 
experimentally determined value o'f -200®F (144K). The experimental 
procedure followed in making this determination is presented in the 
following two chapters. Using classical lamination theory, the trans- 
verse normal stress, C2» predicted in the 90® ply at -200”F is 11.59 ksi 
(79.91 MPa) which corresponds to an in-situ transverse strength of 2.55 
Yt(T). 

The transverse stress, 02* in the 90° layer in the vicinity of a 
transverse crack as predicted by the finite element and modified shear- 
lag analyses are shown in Fig. 13. Finite element results are presented 
at four locations through the thickness of the 90° layer since a through 
the thickness gradient is predicted whereas a single representation is 
given for tne shear-lag prediction which is independent of position 
through the layer thickness. A comparison of the t./o analyses shows 
that the shear-lay analysis predicts the rebuilding of stress in the 90° 
layer to occur in a shorter distance than does the 'Inite element analy- 
sis for the majority of the thickness of the layer. Since a mathematl- 
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cal singularity exists at the tip of the transverse crack in the finite 

element analysis, near the crack tip is greater than the shear-lag 

analysis predictions. Similarly, the shear stress, a- , predicted along 

xy 

the 0°/90° interface is greater in the finite element analysis than in 

the shear-lag analysis (Fig. 14), The shape of the curves for 

both and x are similar for the two analyses and suggest that a 
c xy 

major portion of the stress is recovered in the cracked 90° layer in 

r 

less than 0,05 in. (^ = 5), but that a much greater distance is required 
for ©2 to attain 02 ^^^ ^ major portion of the layer thickness. 

The through the thickness gradients in Oo predicted with the finite 

L. 

element analysis are shown in Fig, 15 at several distances from the 
transverse crack. The gradient in is seen to be most pronounced near 
the transverse crack, decreasing In magnitude with distance from the 
crack plane. As a result of this gradient, material near the 0°/90° 
interface is at the highest level of and therefore is of primary 
interest for purposes of predicting the location of the next transverse 
crack . 

Since the predicted distance at which 02 attains its uncracked 

value, 02 *^^^, in the finite element analysis is dependent upon the 

location through the 90° layer thickness, an assumption m'st be made to 

determine the equilibrium transverse crack spacing. If a crack is 

u 1 1 

assumed to occur when first reaches 02 at any point in the 90° 
layer (disregarding the singularity, i.e., z/T < 0.875 for the mesh 
used) the equilibrium crack spacing predicted by the finite element 
analysis is 0.09 in. (-^ = 9), corresponding to a crack density of 11.1 
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Fig, 15. cr^ Stress Gradients in 90‘' Layer of [02^^^2^s 
Laminate in the Vicinity of a Transverse Crack 
(-200‘’F). 
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cracks/in. (4.4 cracks/cm). If, however, a crack is assume'l to occur 
when Og throughout the entire thickness of the 90° layer reaches 
the predicted crack spacing becomes 0.13 in. = 13) 
corresponding to a crack density of 7.7 cracks/in. (3.0 cracks/cm). 

This represents a 44% difference in predicted crack spacing between 
these two limiting cases. 

In the modified shear-lag analysis, approaches asymptotic- 

ally and therefore a limiting percentage of the in situ strength must be 
chosen to predict a finite crack spacing. Assuming that a transverse 
crack occurs when ^2 reaches 99% of 02^^^» predicted equilibrium 
crack spacing is 0.056 in, = 5.6) corresponding to 17.9 cracks/in. 
(7.0 cracks/cm). If a value of 99.9% of 02 *^^^ is used, the predicted 
crack spacing is increased 50% to 0,084 in. - 8.4), a crack density 
of 11.9 cracks/in. (4.7 cracks/cm). A good correlation exists between 
the tv;o analyses v/hen assuming crack formation occurs when 02 first 
reaches ^ 2 *^^^ in the 90° layer of the finite element analysis 
(disregarding the singularity) and using 99.9% of ^ 2 ^^^ if’ modified 
shear-lag analysis (9.0 y and 8.4 ^ crack spacing, respectively). At 
the temperature of initial transverse cracking (-200°F), no further 
transverse cracking will occur beyond the equilibrium crack spacing in 
either analysis since a, is less than throughout the 90’ layer. 

Knowing the initial equilibrium crack spacing, the next task is to 
determine the additional thermal load needed to produce further trans- 
verse cracking. The location of the next cracks are assumed to be at 
the midplane between existing cracks since 02 Increases with distance 
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from each transverse crack and is greatest at those planes. Thus the 
task of determining the thermal load required to produce additional 
transverse cracking becomes one of determining the temperature at 
which Op at the midplane between existing cracks reaches the in situ 
transverse tensile strength, 

The Op stress in the 90" layer of the [02‘^^^2^s ^^minate at the 
mioplane between two existing cracks (y=b) is shov/n in Fig. 16 for three 
crack spacings. As might be expected, the through the thickness grad- 
ient predicted by the finite element analysis increases In magnitude as 
the crack spacing decreases. At the initial equilibrium crack spac- 
ing (-^ = 9), the laminate midplane (y *• C.O) i s at a 347, lower value 
of o^ than the 0'’/90" interface (y ^ 1.0). Assuming that failure occurs 
when 02 first reaches o^^^^ 'it any point in the 90" layer, a large 
percentage of the layer thickness is well below at the mioplane 

C 

between cracks (y=B), The shear-lag analysis predictions are seen to 
correlate well with the finite element analysis predictions near the 
0®/90" Interface. 

For predicting additional transverse cracking under increased 
thermal loading, an initial equilibrium crack spacing of 0.09 in. (0.23 
cm), or 2B/T = is assumed to be present in the [ 02 / 902 ]^ laminate at 
-2C0°F. In the finite element analysis, additional cracking is assumed 
to occur at the temperature at which O 2 first attains 02 ^^^ at a point 
in the midplane between existing cracks. Because of the temperature 
dependent niaterial properties, an iterative procedure is required to 
detenmine this temperature. For the laminate, this tempera- 
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tore is predicted as -235°F (125K). The crack spacing is thus rediiced 
to one-half the initial equilibrium crack spacing, or 0.045 In, 

(— Y- = 4.5, or 0.11 cm) which corresponds to a crack density of 22.2 
cracks/in. (8.7 cracks/cm). In the modified shear-lag analysis, Eqn. 
(3.30) along with a classical lamination theory analysis are used to 
predict the temperature at which attains © 2 ^*^^ the midplane 
between existing cracks and additional cracking occurs. This tempera- 
ture is predicted as -225®F (130K), 10°F higher than with the finite 
element analysis. Neither analysis predicts any further transverse 
cracking prior to -250“F, the lowest temperature used in the experimen- 
tal portion of this investigation. 

3.4 Influence of Transverse Cracking on Laminate GTE 

A generaliaed plane strain finite element analysis is used to 
determine the influence of transverse cracks on tiie CTE for the ten 
laminates listed at trie beginning of the chapter, C.racks af'e assumed to 
be present at ever, intervals and the resulting modeled region and bound- 
ary conditions are identical to those used in deterr’ning the initial 
equilibrium crack spacing (Fig. 7). To determine the laminate CTE, the 
modeled region is subjected to a uniform temperature change (taken to 
be AT -• - l^F) at room temperature, Tne unknown, uniform displacement, 
V*, of the line y=B is solved for and the laminate CTE, is given by 

“y " 'ea t 


(3.31) 
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where B Is the length of the modeled region (one-half of the cracK 
spacing). The 272 elcirient, 306 node finite element mesh shown In Fig. 

10 Is again used. The transverse crack density was varied by changing 
the y scale factor of the mesh. Each laminate was analyzed for several 
transverse crack densities ranging from 2 cracks/in. (0.79 cracks/cm) to 
100 cracks/in. (39.4 cracks/cm). In addition, uncracked laminates were 
analyzed and the result compared to the classical lamination theory 
solution. 

The variation of laminate CTE with crack density in cross-ply 
laminates is shown In Fig. 17. These same results are presented In the 
form cf percent retention curves in Fig. 18. The curves for the quasi- 
isotropic laminates (Fig. 19) represent both actual changes and percent 
retention since the three quasl-isotrcpic laminates have the same CTE In 
the uncracked state. As expected, the finite elonent results show 
excellent agreement with classical lamination theory for all ten lami- 
nates In the uncracked state. 

The cross-ply laminates exhibit a rather sharp drop in CTE at low 
crack densities (less than 25 cracks/in. (19.7 cracks/cm). The percent 
retention curves of Fig. 18 predict a range of 35-76 percent reduction 
at a density of 25 crccks/ln. (9.3 cracks/cm) and a range of 50-83 
percent reduction at 50 cracks/in. (19.7 cracks/cm). Crack densities as 
high as 78 cracks/in. (30.7 cracks/cm) have been observed in quasi- 
Isotropic T300/5208 graphite-epoxy laminates [8] and as high as 100 
cracks/in. (39.4 cracks/cm) In graphite-polyimide [4]. The percent 
retention curves also show that the change in CTE is largest in the 
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Fig. 17, CTE as a Function of Crack Density in 
Cross-Ply Laminates. 
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laminates with the highest percentage of 90° layers ([ 0/90335 and 
[ 903 /Ojj). For the three pairs of Icmlnates with the same percentage of 
90° layer, the laminate with the 90° layers at the exterior ([90p/0,,,l5) 
exhibit a slightly greater reduction In CTE than with the 90° layers 
together In the interior ([Oj^/OOi^]^) . This result is believed to be due 
to the additional freedom of the [OO^^/O^jj laminate to experience strain 
In the z direction (Fig. 8 ) in the vicinity of a transverse crack. 

The results for the quasl-isotropic laminates (Fig, 19) show that 
the CTE decreases with increasing crack density, b-jt at a nwch lower 
rate and to a lessor extent than the cross-ply laminates of Fig. 10. 
These differences are largely the result of the lower percentage of 90° 
layers In the laminate. However, as Indicated In Fig. 19, the CTE 
retention Is not controlled only by the percentage of 90° layers. It is 
also a function of the stacking sequence of the laminate. The laminates 
with the outer 90° layers or the two 90' s adjacent at the midplane have 
only a 57 percent retention at 100 cracks/in. (39,4 cracks/cm) whereas 
the laminate with 90's interspersed between +45° layers has a 65 percent 
retention at 100 cracks/in. This difference is due to the higher con- 
straint 1n the thinner interior 90 layer. As in the cross-ply lami- 
nates, placing the 90° layers at the exterior of the laminate 
([90/+45/0]j and [90/-45/0/4535) results in a slightly greater reduction 
in CTE than when together at the interior ([0/t45/90j5). These results 
also clearly show that a lamination type theory which is independent of 
stacking sequence effects is not sufficient for analysis of this 
problem. 


Chapter 4 


EXPERIMENTAL PROGRAM 


4.1 Material Specificatio n 

Five panels viere used in this study, all from the same batch of 
prepreg material. The panels were fabricated using Narmco T300/5208 
graphite-epoxy prepreg tape and cured according to the manufacturer's 
recommendations. All five panels were examined for flaws using a 
C-scan ultrasonic detector and were found to be of good quality. 

Fiber volume determinations were made on all panels. The defini- 
tion of a percent fiber volume used was: 


Using the relation = Wj: + where = weight of composite sample, 
equation (4.1) may oe expressed as 


V 



where 


= weight of fiber 
= weight of matrix 
p^ = density of fiber 
p = density of matrix 


m 


V 
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For T300/5208 the following values were used: 


= 0.0636 Ib/in-^ 

= 0.0459 Ib/in^ 


(1.76 g/cm'’) 
(1.27 g/cm^) 


Therefore, equation (4.2) becoines 


V 


100 


1.386 (jj^ - 1) + 1 


(4.3) 


Thus a fiber volume determi nation requires weighing the conposite sample 
(W^), removing the matrix material using nitric acid digestion, and 
weighing the fiber (H^). Fiber volumes for the five panels were deter- 
mined to be within sixty-seven and seventy-three percent (Table i). 


4,2 Specimen Preparation 

A total of twenty-four specimen': were constructed, six from each 
cross-ply panel and three from each ciasi-isotropic panel. All speci- 
mens Mere cut to 2.50 + .01 in. (6.35 + .025 cm) squares. To enhance 
microscopic observation, two adjacent edges of each specimen were 
polished in the following manner. Saw marks were removed and the 
surface made planar using 600 grit paper. After cleaning, the specimen 
was polished on a inetallographic wheel using the final abrasive, a 0.3 
micron diamond paste solution. To avoid rounding the edges, the 
specimen was supported against a piece of aluminum angle while 
polishing. The supported specimen was rotated around the wheel in the 
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Table 1. Fiber Volumes of Test Panels 


Laminate 
Configurations 
in Panel 

Specimen 1 

Specimen 2 

Average 

[O/SOjJj 

[SO/Oj]^ 

68.7 

68.8 

68.75 

[ 02 / 902 ]^ 

[902/02], 

69.9 

71.9 

70.9 

[ 03 / 90 ]^ 

[ 903 / 0 ]^ 

72.6 

71,8 

72.2 

[0/±45/90]3 

[90/+45/0]5 

67.8 

68.7 

68.25 

[0/45/90/-45]3 

[90/-45/0/45]3 

72.4 

72.0 

72.2 
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opposite direction to its motion to avoid a directional bias in the 
polishing. The procedure required about one-half hour for each surface. 

To identify the location of transverse cracks during microscopic 
observation, ttiin tape markings were placed in 0.2 in. (0.5 cm) 
intervals on the specimen surface along the two polished edges (Fig. 

20) . Such markings were visabic as white humps on the specimen edge 
when viewed through the microscope. 

Prior to thermal loading all specimens were dried for at least one 
v'eek (168 hours) in a Hotpack vacuum oven at 30 in. Hg vacuum and 129°F 
(327K) and were stored in the chamber prior to testing and betv;een 
thermal cycles. 

4.3 Equipf.ient 

The development of transverse cracking during thermal loading was 
monitored through a 3 x 6 in. (7.6 x 15.2 cm) window on the side of an 
Applied Test Systems liquid nitrogen cooled environmental chamber (Fig. 

21) , An analog temperature controller used in conjunction with a shut 
down value on the liquid nitrogen supply line provided control of the 
temperature and cooling rate of the chamber. 7ne terr^perature of the 
chamber was monitored with a thermocouple bonded to the surface of a 
control specimen with silicone rubber and suspended directly above the 
test specimen. 

Specimens were held in a rotating holder during testing to enable 
the viewing of either polished edge (Fig. 22). Us-ng tv .'0 nylon cords 
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Fig. 20. Thermal Cycling Specimen 
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Fig. 22. Interior View of Environmental Chamber. 
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extending out through the top of the chamber, the specimen could be 
rotated without opening the chamber door. 

A 75X microscope with a 2.1 in. (5.3 cm) focal length was used to 
detect transverse cracking during a test. In detecting cracking, the 
light source was directed nearly perpendicular to the specimen surface 
but at 3 slight angle to help exaggerate surface relief. After 
completing a thermal cycle, a Unitron microscope was used at higher 
magnifications (100-400X) to confirm cracking located during the test. 
Photography was also possible using a Polaroid camera attachment. 

4 . 4 Initial Thermal Cycle Tests 

Transverse cracking was monitored in 25”F (14K) increments during 
the initial cool down to -250®F (115K) beginning at 75°F (297K). A 
cooling rate of approximately 10“F (5.6K) per minute was used between 
examiration temperatures, at which times the temperature was kept 
approximately constant. Specimens were tested individually due to the 
length cf time required to carefully examine the two polished edges. 
Starting at a corner of the specimen, one edge was slowly traversed with 
the microscope and the location of cracks recorded with respect to the 
tape markings on the specimen. After completel> viewing one edge, the 
specimen was rotated and the procedure repeated on the adjacent edge. 
Cracks were visible as thin lines, usually very straight and perpendic- 
ular to the top and hotter surface (Fig. 23). A complete test required 


from four to six hours 
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Fig. 23. 


Transverse Crack'>ng in a [0^/902]^ Laminate. 
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After the final examination at -250'*F (116K), the specimen v«s 
slowly warmod to room temperature, removed frcm the chamber, and exam- 
ined at higher magnifications (100-400X) to confirm the presence and 
location of cracks detected during the test. 

4 , 5 Thermal Cycling Tests 

Following the initial thermal cycle tests, fifteen specimens (three 
from each uf the five panels) were thermally cycled between -250 and 
+250°F (116 and 394K)„ The initial cool down and return to room temper- 
ature was considered as the first thermal cycle. Subsequent thermal 
cycles consisted of heating from ream temperature to 250“F (394K) , 
cooling to -250‘’F (1J6K), and heating to room temperature all at a rate 
of aeproximately 10*F per m'nute. TV>.e specimens were exposed tc 

a total cf twenty thernial cycles. Transverse cracking was monitort'd 
with a Uhitron microscope (100-4.10X) at tie completion of the first, 
second, third, fourth, fifth, se^'enth, tenth, fifteenth, and twentieth 
cycle. 


Chapter 5 

EXPERIMENTAL RESaiTS AND OHSERVATIONS 

As described in sections 4.4 and 4.5, two types of tests were 
performed on the T300/5208 graphite epoxy specimens. Initial thermal 
cycle tests provided a knowledge of the temperature at which transverse 
cracking first occurred in a laminate and the crack density as a func- 
tion of temperature. Tliermal cycling tests provided a knowledge of 
crack density as a function of the number of cycles. In addition, 
several observations regarding the characteristics of transverse crack- 
ing were made while performing these tests. This chapter presents a 
discussion of all experimental results and related observations. 

Prior to testing, all specimens were microscopically observed and 
found to be free of transverse cracks. The outer layers of the speci- 
mens cut from the t'vo quasi -Isotropic panels, however, were found to 
have large amounts of matrix damage present in the form of randomly 
directed crazing (Fig. 24). Consequently, an analysis of transverse 
cracking ir. the 90° layers of the [90/+45/0]^ and [90/-45/0/45J5 lami- 
nate configurations was not possible. This damage was believed to have 
occurred either while cutting the specimens or during surface prepara- 
tion. 

Transverse crack densities reported are average densities over the 
length of the specimen (2.5 in. (6.35 cm)). Crack densities in the 
laminate configurations wich two separate 90° layers ([9C/O235. 
[902/0235, 4nd [9O2/035 laminates) are given as averages of the tv/o 
layers. 
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a) [90/-45/0/45]^ Laminate 



Fig. 24. Outer 90° Layer D.amagc Present Pt ior to 
Testing in Quas i -Isotropic Laminctes. 
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5.1 Initial Thermal Cycle Test Results 

The temperatures at which transverse cracking was first observed in 
the cross-ply laminate configurations are shown in Table 2. During the 
initial cool-down from room temperature to -250°F (116K), four of the 
six cross-ply laminates experienced transverse cracking. A total of 13 
of the 36 cross-ply specimen configurations exhibited transverse crack- 
ing. Of these 13 specimen configurations, cracks were observed during 
the test In only 8 specimen configurations; cracks were observed in the 
remaining 5 after the tests were completed and at higher magnifications. 

A comparison of tne three cross-ply laminate configurations with 
interior 90° layers ([O/QO^l^, [02/90,]^, and [0^/90]^) is shown in Fig. 
25. As the thickness of the 90° layer decreased from the [O/QO^jj 
laminate to the [Oj/OO], laminate, the thermal load required to produce 
transverse cracking decreased. This trend follows the classical lamina- 
tion theory prediction (Chapter 3) and is a result of the corresponding 
Increase in the 0° layer thickness which produces an increased con- 
straint on the 90° layer and therefore an increased transverse normal 
stress, 02* at a given thermal load. 

A comparison of the temperatures of initial transverse cracking in 
the three cross-ply laminaces with exterior 90° layers {[902/0]^, 
[902/'02]5, and [90/0335) is shown In Fig. 26. Following classical 
lamination theory analysis, it was expected that transverse cracking 
would first occur in the [50/0335 laminate Followed by the [502/023^ 
laminate and then the [9.13/035 la.minate. Although the [9O/O335 laminate 
was the first to cxperieice transverse cracking, the [90'j/035 was the 




Table 2 


Temperature of Initial Transverse Cracking, ®F (K) 


Laminate 

Configuration 

specimen Number 

1 

Initial 
Cracking 
Temp . 

1 

2 

n 

3 

4 

5 

6 

[ 0/90315 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

[9O/O3I5 

-50 

(228) 

UD 

UD 

NC 

UD 

NC 

-50 

(-228) 

[ 02/90215 

-200 

(144) 

NC 

NC 

UD 

UD 

-200 

(144) 

-200 

(144) 

[9O2/ ^ 2-5 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

[ 03/9015 

NC 

NC 

-175 

(158) 

NC 

NC 

NC 

-175 

(158) 

[ 903/015 

-75 

(214) 

-175 

(158) 

-100 

( 200 ) 

NC 

-225 

(130) 

NC 

-75 

(214) 


NC - No cracking. 

DO - Cracking undetected during test, confirmed after test. 


ORIG!;\'AL F.M ~ I'J 
OF POOR QUhLITV 


69 


oi^lCi^'Al. r;.cz 'S 

OF POOR QUALllV 



Fig. 25. Temperature of Initial Transverse Cracking 

in 90° Layer of [0 /90 ] Laminates, 
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next tc crack and the [902/02]s laminate remained crack- free for all 
tests. All transverse cracking in the [90^/0]5 laminate prior to -225°F 
(130K) occurred in specimens 1, 2, and 3 all of which experienced 
sudden, extensive transverse cracking with noticeable vyarping (Fig, 27) 
at increased thermal load levels (-200, -175, and -225°F, respec- 
tively). Microscopic observation revealed wide transverse cracks often 
accompanied with delaminations, primarily in the layer at the outside of 
the curvature (Fig. 28). Some fiber breakage in the adjacent 0® layer 
was also observed (Fig. 29). 

While three [903/0]g laminate specimens experienced extensive 
transverse cracking and v/arping, the remaining three specimens tested 
were found to be more resistant to cracking and did not experience 
warping. This difference in behavior is believed to be a result of the 
position of the specimens In the panel from which they were cut. As 
shown in Fig. 30, specimens 1, 2, and 3 (those experiencing extensive 
cracking and warping) were cut from a portion of the panel adjacent to 
the edge after trimming off 0.5 in. (1.3 cm) whereas specimens 4, 5, and 
6 (those more resistant to cracking) came from an interior portion of 
the panel. An investigation of the unused portion of the panel showed 
variations in the laminate thickness as far as 1.5 in. (3.8 cm) inward 
from the edge. Therefore, a thickness variation may have been present 
in specimens 1, 2, and 3. Although these specimens appeared to be 
uniform in thickness and symmetric under microscopic observations, it is 
believed that these variations were at least partially responsible for 
the observed behavior. Since the same specimens were used to study both 
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Fig. 27. Warping of [CO^/O]^ Laminate Specimens 
1, 2, and 3 Due to Thermal Loading. 
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Fig. 29. 


Fiber Br•.....^age in Adjacent 0° Layer of 
[ 902 / 0 ]^ Laminate (Specimen 2). 
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the [90j/03j. and laminates, some caution should be used when 

interpreting the [0j/903j laminate results, especially in the cyclic 
thermal tests where only specimens 1, 2 , and 3 (those experiencing 
warping) were tested. Five of the six specimens, however, remained 
uncracked in the [C.,/903^ laminate configuration through the initial 
thermal cycle tests. 

It is also noted that specimens 1, 2, and 3 used to study both the 
C0/90j]j and [90/0^]^ laminates were cut from the edge of the panel 
after trimming off 0.5 in. (1.3 cm). However, no unusual behavior 
occurred from specimens 4, 5, and 6 which were cut from interior 
portions of the panel. All other specimens were cut from interior 
portions of the panels. 

After the initial formation of transverse cracking was observed, 
additional cracking was reported under increasing thermal load in three 
cross-ply laminate conf i^rations as shown in Table 3. At the co't^le- 
tion of the initial cool down to -230°F (1161;), all specimens were 
SiOwly wanned to room temperature and viewed under a microscope at 
higher magni f icstions. Numerous transverse cracks, undetected du-ing 
the test, w^re obs*»rved in both the [ 02 / 902 ]^ and [90/9^]^ laminates. 
Failure to detect cracks in the Cfl2'^^2^s l^'‘”ihate was believed to be 
due to the small ci'ck opening whereas slight rounding of the specimen 
edges during polishing made crack oetection difficult in the thin 
[0/90j]j. laminate outor layers. 

At the completion of the initial cool down, the highest crack 
density was present in the [90-j/0]j. specimens with warping present 


Table 3 


Transverse Crack Density Versus Temperature 
cracks/in. (cracks/cm) 


Laminate 
Conf i g>jr- 
atior. 

Tenperaturg. °F (K) 

Crack 

Density 

After 

Test 

0 0 

70 

•V 0 
0 5 * 
0 

^ t- 

0 -j 
c:. . 

:.L 
r- M 

-^v 

-50 

(228) 

-75 

(214) 

-100 

( 200 ) 

-125 

(186) 

-150 

(172) 

-175 

(158) 

-200 

(144) 

-225 

(130) 

-250 

(115) 

[D/90^J5 

- 

- 

- 

- 


- 

- 

- 

- 

- 

[ 90/0315 

0.0333 
(0.131 ). 

0.0333 

(0.0131) 

0.066? 

(0.0262) 

0.056/ 

(0.0262) 

0.0667 

(0.0262) 

0.0667 

(0.0262) 

0.0667 

(0.0262) 

0.667 

(0.0262) 

0.0667 

(0.0262) 

0.233 

(0.0919) 

[02/902]s 

- 

- 

- 

- 


- 

0.200 

(0.0787) 

0.200 

(0.0787) 

0.533 

( 0 . 210 ) 

1.27 

(0.499) 

[ 9 O 2 /O 2 J 5 

- 

- 

- 

1 

1 . 

- 

- 

- 

- 

- 

Warping 

- 

1 


! 

! 

- 

- 

- 

- 

- 

[ 903 / 0 ] 

No 

Harping 

- 

- 

1 

1 

i 

1 - 

i 

1 

- 

- 

- 

0.0667 

(0.0262) 

0.0667 

(C.0262) 

0.0657 

(0.0262) 

[ 03/0015 

Warping 

- 

- 

- 

i 

- 

0.133 

(0.0525) 

0.133 

(0.0525) 

0.133 

(0.0525) 

0.133 

(0.0525) 

0.133 

(0.0525) 

[ 903/015 

Warping 

- 

0.0667 

(0.0262) 

0.467 

(0.184) 

0.667 

1(0.263) 

0.800 

(0.315) 

1.67 

(0.656) 

2.20 

( 0 . 866 ) 

3.33 

(1.31) 

3.33 

(1.31) 

3.40 

(1.34) 
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(specimens 1, 2, and 3). V(1th the exception of the [C'2/902]5 laminate, 
all other laminates had a final crack density of less than 0.25 
cracks/ln. (O.iO cracks/cm) corresponding to an average of less than one 
crack per specimen. The location of transverse cracks appeared to be 
random throughout the initial cycle tests for all laminates and a uni- 
form spacing was not attained in any specimen at the completion of one 
cycle of thermal loading. 

5,2 Cyclic Thermal Test Results 

After the Initial cool down to -250®F (116K), specimens 1, 2, and 3 
from each of the five panels were thermal cycled 20 times between -250 
and 250®F (il6 and 394K). The effects of the cyclic tests on transverse 
crack densities are shown in Table 4. 

5.2.1 Cross-ply Laminates 

At the completion of 15 thermal cycles, transverse cracking was 
present in all six •.'oss-ply configurations. Crack dansities ranged 
from 8.27 cracks/in. (3.26 cracks/cm) in the [03/90],. laminate to 0,133 
cracks/ln, (0.052 cracks/cm) in the [0/903]^ laminate. A comparison of 
the three cross-ply laminate configurations with interior 90“ layers is 
shown in fig. 31. As the thickness of the 90“ layer decreased from the 
[0/903 ]j. laminate to the [03/90]^ laminate, the transverse crack density 
present at a given .lumber of thern’.al cycles increased. As was the case 
with the initial thermal .‘.yde test results for the same three lami- 
nates, this trend follows the classical lamination theory predictions of 


Table 4 

Transverse Crack Density Versus Number of Thermal Cycles 
cracks/in. (cracks/cm) 


Laminate 

Configuration 

Number of Thermal Cvcles, to -250 (115 to 394 K) 

1 

2 

_ 

4 

5 

7 

10 

15 

20 

[0/90, 3s 

- 

- 

- 

- 

- 

- 

- 

0.133 

(0.0525) 

0.133 

(0.0525) 

[ 9 O/O 3 J, 

0.400 

(0.158) 

1.27 

(0.499) 

1.60 

(0.630) 

2.47 

(0.971) 

2.60 

( 1 . 02 ) 

2.87 

(1.13) 

3.60 

(1.42) 

4.07 

(1.60) 

4.07 

(1.60) 


1.07 

(0.420) 

1.33 

(0.525) 

1 .73 
(0,682) 

1.73 

(0.682) 

1.87 

(0.735) 

2,00 

(0.787) 

2.53 

(0.997) 

3.20 

(1.26) 

5.20 

(1.26) 

[POp/Oj], 

- 

- 

- 

- 

- 

- 

0.133 

(0.0525) 

0.200 

(0.0787) 

0.200 

(0.0787) 

1 — J 

0 

CO 

0 
1 — 1 

1 

0.133 

(0.0525) 

1 .733 
fO.682) 

2.27 

(0.892) 

2.93 

(1.15) 

3.07 

( 1 . 21 ) 

5.73 

(2.26) 

6.00 

(2.36) 

8.27 

(3.25) 

8.27 

(3.25) 

[ODp/ri], 

3.33 

(1.31) 

4.33 

(1.71) 

4.33 

(1.71) 

4.33 

(1.71) 

4.33 

(1.71) 

4.60 

(1.81) 

4.67 

(1.84) 

4.67 

(1.84) 

4.73 

( 1 . 86 ) 

[0/+45/90], 

- 

- 

- 

- 

- 

- 

0.133 

(0.0525) 

0.133 

(0.0525) 

0.267 

(0.105) 

[DO/fAS/Ol, 

N/A 

N/A 

N/.A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

[0/45/9D/-4535 

- 

- 

- 

- 

- 

- 

- 

- 

- 

[90/-46/0/45j^ 

f!/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 


N/A - Information Net Available Due to Damage Present Prior to Testing. 
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Fig. 31. Cyclic Thennal Test. Results for ^0 /90 ] 
Cross-Ply Laminates. 
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an Increasing tiansverse normal stress, o^, in the 90° layer with 
decreasing 90° layer thickness at a given thermal lead. The differences 
In crack density Increased v;lth the number of thermal cycles up to 15 
cycles at which point no additional cracking occurred. Throughout the 
thermal cycling, the location of transverse cracks still appeared to be 
quite random and without regiilar spacing. 

A comnarlson of the three cross-ply laminates with exterior 90° 
layers is shown in Fig. 32. The [90^/0]^ laminate, which experienced 
extensive transverse cracking and warping In the initial cycle, showed a 
42% Increase In crack density after 20 cycles. This final crack den- 
sity, however, was only 16% greater than the final •"'anslty of the 
[ 9 O/O 3 IJ laminate, which was nearly crack free after the Initial thermal 
cycle tests. Although predicted to have a higher Oj, In the 90° layer at 
a given thermal load than the [ 9 O 2 /O],. laminate, the laminate 

remained more resistant to transverse cracking. This may have been a 
result of the warping present in the [ 903 / 0 ]^. specimens as discussed In 
the previous section. As was the case wih the internal 90° layer 
cross-ply laminates, c 'ack spacings were .till irregular after 20 ther- 
mal cycles, suggesting that further transverse cracking may have 
occurred with additional thermal cycling. 

An interesting comparison may be made between the crack densities 
in the [ 03 / 90 j^_ and [ 90 / 03 ]^ laminates (Fin, 33). In the uncracked 
state, lamination theory predicts to be identical In the two lami- 
nates at a given ther’al load since tney have the same percentage of 
0 °and 90° layers. !t might therefore be expected that the crack densi- 
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(116 to 394 K) 


Fig. 32. Cyclic Thern;al Test Resiilts foe [90 /O J 
Crcss-Ply Laminates. ' ' 
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Fig. 33. Comparison of [0,/80]_ and [90/C..]. Laminate 
Cyclic Thermal Test Results. 
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ties would be approximately equal after cyclic loading. After 5 thermal 
cycles, however, the crack density in the [03/90]^ laminate was approxi- 
mately twice that in the [90/03]^. indicating an influence of the lami- 
nate stacking sequence on transverse cracking. This difference in crack 
densities may possibly be understood by considering an equal number of 
critical flaws producing transverse cracking to be present in the 90 ° 
layers cf the two laminates per unit volume. A flaw in the [03/90]^ 
laminate would produce a transverse crack extending completely across 
the 90 ° layer whereas a transverse crack produced in the [90/03]^ lami- 
nate would be present in only one of the two exterior 90 ° layers separa- 
ted by the 0 ° layer. The resulting average crack density present in the 
[03/90]^ laminate would be twice that present in the [9Q/O3]. laminate, 
a fact which agrees with the cyclic thermal test results after five 
thermal cycles. 

Similar reasoning suggests that the transverse crack density pre- 
sent in the [03/902]^ laminate should be twice that present in the 
[90,702]^ laminate. Inrough 20 thermal cycles, how'ever, the crack 
density present in the [ 02 / 90 ^]^ laminate was over an order of magnitude 
greater than in the [902/0215 laminate. Considering that such a trend 
developed in the [03/90]^ and [90/0215 laminates only after 5 thermal 
cycles and chat the 90 " layers are at a higher stale of stress, at a 
given thermal load, than in the [O2/9O2I5 and [9O2/O2I5 laminates, it is 
suggested that this trend may have developed if more thermal cycles were 


considered. 
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A similar comparison of the [ 0 /? 03 j 5 and [gOj/O]^. laminates is 
difficult. since the [90^/0]^ laminate specimens used ware those which 
experienced extensive transverse cracking and warping in the initial 
thermal cycle tests. It is believed that the trend of twice the crack 
density present in the laminate with an interior 90° layer might only be 
seen after further thermal cycling, since the state of stress in the 90° 
layers at a given thermal load is the lowest of all the cross-ply lami- 
nates studied. 

5.2,2 Quasi-Isotropic Laminates 

As previously discussed, analysis of transverse cracking in the 
outer layers of the quasi -isotropic panels ([90/+45/0]^ and [90/- 
45/0/45j laminates) was not possible die to extensive matrix damage 
present prior to testing. The [C/+45/90]^ laminate experienced limited 
transverse cracking beginning after ten cycles. No cracking was repor- 
ted in the [0/45/90/-45]j laminate or in any of the +45° or -45° layers 
of any laminate. While the quasi-isotrcpic laminates were found to be 
generally resistant to transverse cracking, a slight increase in 
resistance was exhibiten when the stacking sequence was changed from 
[0/+45/90]j to [0/45/90/-45]g resulting in a thinner 90° layer. 

5.3 Correlation with Ai-jalytical Studies 

5.3.1 Determination of In Situ Transverse Strength 

Experimental results for the initial temperature of transverse 
cracking along with a classical lamination theory analysis (Section 
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3.2.1) may be used to determine the in situ transverse strength, 
u 1 

02 of the laminates investigated. In making these determinations, 

the experimental ly determined temperature of initial transverse cracking 

in a laminate is input into a classical lamination theory analysis and 

the resulting transverse stress in the 90° layer 
u 1 ^ 

(©2 ) determined. Results for the six cross-ply laminates are shown 

in Fig. 3^, All cross-ply laminates exhibit an in situ transverse 
strength, at least 1.9 times greater than the uniaxial 

transverse strength Yj.(T) . Some specimens frem all six laminates 
remained uncracked at -250°r (116K), indicating that Og'*^^ n:ay be 
somewhat greater than the mvaximum values designated by dashed lines in 
Fig. 34. 

II 1 t 

A 90° layer thickness effect on 02 niay be seen in the [0^/90^]^ 
tross-ply laminates. As the thickness of the 90° layer decreases, 

02 ^^^ increases, Trils result is in agreement with ref. [34] and may be 
a result of the corresponding increase in the constraint Imposed by an 
increasing thickness of adjacent 0° layers. A similar thickness effect 
was not evident in the [90^/0^]^ cross-ply laminates. 

Since all four quasi -Isotropic laminates remained iincracked 
throughout the initial cool down to -250°F, only a lower limit for 
may bo determined. At -250°F, - 2.66 in the 90° 

C C '■ 

layers of the four quasi-isotropic laminates. 

In determining in situ strength using experimentr'l results and a 
classical lamination theory analysis, it is noted that finite length 
specimens with free edges v^ere used in the expcr1tr,ental study whereas 
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Fig, 34. In Situ Transverse Strength of 90'’ Layers in Cross-Ply Laminates. 
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classical lamination thaory assumes infinite length laminates Nagarkar 
and Herakovich [32] studied fi'^'te width composite laminates and found 
the thermal st'.'ts> o^, in the 90“ layer of a [0/90]c T300/5208 laminate 
to be 26 percent greater at the free edge than in the interior portion 
of the laminate. For a [90/0]^ laminate, edge was If* 

greater than in the interior. Therefore, the temperature of initial 
transverse cracking may correspond to a crack forming at the free edge 
of the specimen whereas the stress predicted with classical lamination 
theory is the lower value present in the interior of the rpecimen. 

Thus, the in situ transverse strength, even greater than 

predicted. 

5.3.2 Crack Spacing in [02/902J5 

A comparison of the crack spacing observed exparimontal ly (section 
5.1) ano predicted analytically (section 3.3.3) for the [82/'*^2^s 
laminate during the initial cool down to -250“F (116K) is shown in Fig. 
35. The uniform crack spacings predicted with a finite clement analysis 
(0.09 in. at -200°F and 0.045 in. at -250°F) are great contrast to 
the sparse, randomly located transverse cracking observed 
experimentally. As shown in Fig. 36, transverse crack spacing was still 
non-uniform throughout the specimen after 20 thermal cycles. However, 
some localized regions had developed approximately uniform crack 
spacings which approached the predicted spacing. It is believed that 
the predicted uniform crack spacing may have developed throughout the 
length of tlie specimen with additional thermal cycling. 
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Fig. 35. Transverse Cracking Cnjring Therir^al Loading In [02'^^®2^s 
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5.4 Observations 

i)ur1rig both the Initial thermal cycle tests and cyclic thermal 
tests, all transverse cracks were found to be fully propagated across 
the thickness of the 90“ layer upon first fiotection. The location of 
the Initial transverse cracks appeared to be random, favoring neither 
the end nor the midpoint of the specimen. Transverse cracks tended to 
avoid resin rich regions, and as shown In Fig. 37, generally changed 
direction to follow the edge of the region rather than pass through it. 

As may be seen in F1g. 38, transverse cracking involved fiber 
splitting in addition to matrl/c cracking. Generally, fiber splitting 
provided for a straight crack path, seldom meandering from Its previous 
direction. Fiber splitting was observed In all laminates exhibiting 
transverse cracking. 

Transverse cracks were seen to vary in width across the 90“ layer, 
tapering when approaching an adjacent layer as shown In l-ig. 38. Sjch 
tapering was present In all laminates exhibiting cracking and Indicates 
the constraint placed on the cracked 90“ layer by the adjacent layer. 

The presence of delarnl nations acconpanying transverse cracks was 
observed in the [903/0]^ laminate specimens experiencing warping as 
previously discussed (Fig. 28). Such delarnl nations generally followed 
the interface between adjacent layers, but as shown in Fig. 29, some- 
times fiber breakage occurred in the adjacent 0“ layer, 'lel ami nations 
were sometimes found to grow in length with progressive thermal cycling. 

Some transverse cracks In the cross-ply laminates with Interior 90“ 
layers were seen to branch into tw'p distinct crack's as shown in 
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Fig. 37. Transvef-;e Cracking Around a Resin Rich 
Region in a [02/902]^ Laminate. 
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Fig. 39. Generally, the point at which the branching occurred was only 
a few fiber diameters from the adjacent 0° layer. !!o crack branching 
was seen in any of the remaining laminates Investigated. 
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b) [0..,/90]^ Laminate 
Fig, 39. Transverse Crack Branching. 


Chapter 6 
CON’CLUSIONS 

This study has Investigated the characteristics of thermally- 
induced transverse cracks, both experimentally and analytical ly, in 
T300/5208 graphite-epoxy laminates. Of particular Interest was the 
determination of the temperature at which transverse cracks initially 
form, the state of stress in the vicinity of a transverse crack, the 
Initial equilibrium transverse crack spacing, the formation of addi- 
tional transverse cracking under increased thermal load and under ther- 
mal cycling, and the Influence of transverse cracking on the laminate 
CTE. The following conclusions may be drawn as a result of this study: 

1. The T300/5208 laminates studied in this investigation are free 
of cracks after curing. 

2. Four of the six cross-ply laminates experience transverse 
cracking during the Initial cod do»^'n to -250'’F (116K) whereas all 
quasl-lsotropic laminates remain crack free. 

3. Increas1n>^ the percentage of the 90° layer In the L0|^/90 jjj]^ 
laminate results in an increased thermal load required to produce trans 
verse cracking. 

4. The in situ transverse strength of the 90° layer is more than 
1.9 times greater than the transverse strength of the unioi rectional 90 
material . 

5. The initicl equi 1 i briom transverse crack spacing in the 
[ 02 / 902 ]$ laminate is quite random and much greater tnan the uniform 
crack spacing predicted analytically. 
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6 . A through the thickness gradient in the transverse tensile 
stress is present in the cracked 90® layer, decreasing in magnitude with 
distance from the transverse crack and resulting in a higher stress near 
the adjacent constraining layer than in the interior portion of the 90® 
layer. 

7. The initial equilibrium transverse crack spacing predicted in 
the [ 02/90215 laminate with a generalized plane strain finite element 
analysis corresponds well to that predicted with a modified shear-lag 
analysis. 

8 . All six cross-ply laminates as well as the [0/+45/90J5 quasi- 
isotropic laminate exhibit transverse cracking following 20 thermal 
cycles between 250 and -250®F (116 end 394K). 

9. A uniform transverse crack spacing is not present in any of the 
laminates studied following 20 thermal cycles, but the uniformity of 
crack spacing increases with increasing thermal cycles. 

10. The cross-ply laminates exhibit a rather sharp drop in lami- 
nate C 'E at low crack densities (less than 50 cracks/1 n. (19.7 
cracks/ cm)). 

11. The quasi-isotropic laminates exhibit a decrease in laminate 
CTE with increasing crack density, but at a much lower rate and to a 
lesser extent than the cross-ply laminates. 

12. While transverse cracking was able to be detected in a major- 
ity of the laminates with a 75 X microscope, higher magnification is 
needed to confirm the presence of cracks in all laminates studied. 
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Table A.1 f'latorial Property Coefficients for T300/5208 Graphite-Epoxy 

Second Order Least Squares Fit Coefficients to b? used in the Equation: 


Property (T) = + C^T + (T in ^F) 


Property 

^'0 

r* 

^1 

C2 

Reference 

E, (Psi) 

2.13 X 10^ 

-5.S7 X 10^ 

2.32 X 1q2 

25 

i^2.F2 C'sO 

1.63 X 10^ 

-4.10 X 1q2 

2.68 X 10° 

26,27 


1.06 X IC*’ 

-6.00 X 10^ 

-1.22 X 10° 

26 

*^23 

4.87 X 10^ 


1 

3-2 

'■’12 ’■‘'13 

3.43 X 10"’ 

-5.65 X 10"^ 

1.22 X 10"’^ 

26 

'’23 

4.90 X 10"’ 

“ 

- 

32 

xj (Psi) 

2.01 X in^ 

1.05 X 10^ 

1.36 X 10"^ 

26 

yy (Psi) 

4.04 X 10^ 

-1 .73 X 10"’ 

1.17 X 10*2 

26 

U/°F} 

-1.24 X 10"^ 

1.34 X 10"” 

5.45 X 10*’2 

23 

<xp (eZ-Op) 

. 

1.18 X 10"'^ 

i 

1.35 X 10"® 

— j 

2.07 X 10"” 

29,30,31 
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APPENDIX B 


Constrained Displacement Finite Element Solutions 


The process of constraining displacements requires a technique of 
systematically modifying the global stiffness matrix and the global 
force vector of the finite element solution. Branca [40] and Crane and 
Adams [41] have illustrated the modifications required. Hov.-ever, it 
does not appear that a rigorous development of the procedure has been 
presented. Of particular Interest Is the question of a unique solution 
while permiting an arbitrary choice of the constrained node used in the 
solution. This appendix addresses this question in a rigorous manner by 
consioering the variational forntjlatlon of the anisotropic, therrno- 
elastic, generalized plane strain finite element analysis. 

The principle of virtual ivork may be v/ritten 


6(U + y) = 0 


(B.i) 


where 


U = stored strain energy 
V H potential of applied loads 


Equation (B.I) can be written 



N 

(o*c)dR - E 


where 


e = mechanical strain 


Q. = applied nodal forces 


q^ = nodal displacements 
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n 1 

< • » w 


QOAI.i 


; > 
T'i 


N = number of nodes 
R H region of interest. 


Assuming superposition of thermal and mechanical strains, the strain- 
displacement relations in the material coordinate system for a tempera 
ture change AT are: 


Where a^. 
of individual 



a , a are the nonzero coefficients of thermal expansion 
xy’ z 

layers. Using elemental displacements of tne form 


N 

u = T. 4’^(y,z)U. + ^ .X 

i=l 

N 

v= ^ 4<^(y,z)V. 


N 

w = E <i>- (y,z)w. 
i=l ’ ’ 


where 

4 !j(y,z) ^element shape functions 
H nodal displacements 
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= total uniform x strain 
N = number of nodes per element 
The strain displacement relations become 


fe) = CBKq} + {Eq) 


(B.4) 


where 


[Bj H matrix of element shape function derivatives. 



and it is assumed that lamina are in x-y planes. The constitutive 
relation in global coordinates can be written 


{c> = [C] {e} 


(8.5) 


combining (B.4) and (B.5), 

(a) = [C]([3]{q} {£^}) (B.6) 

Substituting Kqns. (R.4) and (B.6) into Eqn. (8.2) yields 


!15 


OF peer; OuALii'/ 



(Cl3](q}+ i {e^))c)R - {q {Q )1 = 0 


(B.7) 


At this point, we consider a one element example (Fig. Bl). For 
this example, the quantity [B]{q) may be written 


CBllq) = 


0 

0 

0 

0 


0 


0 

0 

0 


0 


0 

0 


" 2 ,y 

0 '4^ , 0 


^,y 0 




0 


^'i.y ' 


0 

0 


^ 2 ,z 


"2,z '^2.y 
0 0 (!>. 


0 0 

‘i.y ° 

0 ■«S.. 


*2,y “ 


0 


3,z 

y 


0 




0 

0 


\ 


u 


1 


2 1 


^3 


W3 I 


(B. 8 ) 


Next, the v displacements of nodes 2 and 3 are const rai 'ed (V '2 = V 3 ). 
Under this constraint, either V 2 or may be expressed in terms of the 
other displacement, reducing the n^ oer of independent displacements by 
one. A combination of terms is then possible in [B] and {Q}. 

Selecting V 2 as the "master" displacement representing all constrained 
displacements yields 
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Fig. B.l One Element Finite Element Example. 



117 



With the constraint = V 3 imposed, we proceed to take the first vana- 
tion of the nocilfied equation (7) 

[/* CB*]^[C][3*] dR]{q*) + [B*1 [f]‘{c^idR - {0* > = 0 (B.in) 

Jr ^ 





The resulting set of equations has one less equation than before inipos- 
ing the constraint. Assembling Eqn. (R.IO) in a mere conventional form 


and {F*) « {Q*) - / {e„ }dR (B.ll) 

Jr 

A comparison of the modified quantities [K*] and {F*} with the original 
quantit ies [Kj and {F } suggests a series of systematic ntodi 'l cations, 
convenient to the existing computer code, to acco(r.T,odate constrained 
displacenier^ts. Again, consider the one element example she-.-fn in Fig, 
ij.l with the V'2 and V3 displacements constrained. Cefore applying the 
constraint, the governing finite element, equations are of the form 


where 


[K*]{qM “ {F*} 

r - 

[K*] = / [B*] rc3CB*]0R 


^11 ^12 ^13 '^•15 '^17 *^18 h.9 

Koo f -23 *^24 ^25 *^26 *^27 ®^28 ^‘29 

■^13 ^23 ‘'33 *^34 ^35 ^36 '^37 *^38 *^39 

*^14 *^24 *^34 ^'44 *^45 *^46 *^'47 *^48 *^'49 

^15 4s *^35 *^'45 *^55 ^56 *^57 ^58 ^9 

*^16 ^26 *^36 ^46 *''56 ^66 *^67 *^68 “^69 

*^17 ^27 '^37 ^'47 ^57 *^67 h ? *^78 *^79 

“^la ^28 *^38 'M8 ''58 '^73 ''88 '''89 


Ki9 K.^g i<39 K^g Kg, 3 Kjg Kgg Kgg 



(B.l?) 
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When displacements V2 and V, are constrained, the fifth and eighth 
columns of [K] ore added together term by term and the sum entered as 
either the fifth or eighth column, depending on which displacement is 
chosen as the master displacement. The other column is modified by 
placing zeros in all entries. Mext, the fifth and eighth rows are added 
term by term and the sum placed as the fifth row. The eighth row is 
modified by placing a zero in all entries except a one on the main 
diagonal. Similarly in the force vector, the sum of the fifth and 
eighth entry is placed as the fifth entry and a zero is placed as the 
eighth entry. The final form appears as 



K ,2 

^13 

'■'14 

(IC,5»XT8) 

*16 

*17 

G 

*1S 

i 

l.\ 

/ 

/ 

''IZ 

•^22 

^23 

*^24 

(IC25«2g) 

*26 

*27 

C 

*29 


’'l 



*^23 

»^33 

*"'34 

^'’■35*'''38^ 

*36 

*37 

0 

*39 


i“l 1 

1 




^44 


*46 

*47 

C 

*49 


r^’f 

I 

(K, j'HCjg) 




/S5*’'S8 , 

^*66**63^ 

^*57**78^ 

0 

^*59**S9^ 


[h' 
\ ( 

7- / 

^6 

*^26 


^'46 


*66 

*67 

0 

*69 


r\ 

1 


*'27 

h? 

^47 

'h7*ha^ 

*67 

*77 

0 

*79 


l*>3 

0 

0 

0 

0 

0 

0 

0 

\ 

0 






K.39 

^9 

(Kjg^Kgg) 

*69 

*79 

0 

*S9 


i ‘'i 
\ 1 

' \ 

or 



CK*]lu*} 

= iFM 





(B 

.13) 


Displacement V-^ may have been chosen as the master displacement in which 
case the combining of terms would have occurred in row and column 
eight. Tne solutions obtained when using different Mster displacements 
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are Identical, making the choice arblttary. Since the modified stiff- 
ness matrl);, [K*] In Eqn. (B.13) remains symmetric, the solution 
procedure Is unaffected. Although not evident In the one element exam- 
ple considered, the modifications made 1n the stiffness matrix [K] may 
produce a greater bandwidth In [K*], requiring that the bandwidth be 
recomputed after such modi float ions. For multiple element cases, the 
modifications may be irade on the assembled global stiffness matrix 
rather than modify'ing the affected element matrices separately. 
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APPENDIX C 

Composite Laminate Finite Element-Two-Dimens1ona1-Cor.strfiinecl 
Displacements (CLFE2DC) 

Input Data Sequence 


Cards 1 & 2 
cTiunn 
~TM- 


(20A4) 

Title 


Title Cards 
Contents 


Card 


3 (1615) 
Tolumn 


Master Control 
Contents 


Card 


Card 


Card 


TT 

NEM 

fJumber of Elements 

6-*10 

NODS 

Number of Node Points 

11-16 

NPE 

Number of Nodes Per Element 


= 4 

4 Node Linear Elements 


= 8 

8 Node Serendipity Elements 


» 9 

9 ftode Laorangian EleiT:ent;s 

16-20 

NAflG 

Number of Different Angles 

21-25 

NSOr 

Number of ^eclfled Itodal Dlsplccements 

26-30 

NSBF 

Number of Specified Nodal Forces 

31-35 

NCON 

Number of Constrained Displaceflierits 

35-40 

NPLOT 

Plot Option 


= 0 

none 


= 1 

undeformed plot only 


^ 2 

deformed plot or ly 


- 3 

both undeformad and deformed 

41-45 

NTl 

Dump file for displacements 


= 0 

no dump 


^ 0 

dumps displacements to unit NT! 

46-50 

NT2 

Dump file for strains, stresses, and 
strain energies 


= 0 

no dump 


/ 0 

dumps to unit NT2 

51-55 

NCHECK 

Data Check Option (NE, 0 - Data Check Only) 

4 (8D10.5) 

Scale Factors 

tolumn 


Contents 

l-lo 

sew 

Y-Nodal Coordinate Scale Factor 

n-20 

SQM 

Z-Nodal Coordinate Scale Factor 

5 (1615) 

Printer Control Card 

Column 


Contents 


KEY(l) 

Key for Printing Flement Data 

6-10 

KEY (2) 

Key for Printing Nodal Data 
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11-15 KEY (3) Key for Printing Specified 

Displacements 

16-20 KEY (4) Key for Printing Specified Nodal Forces 

21-25 KEY (5) Key for Printing Displacements 

26-30 KEY (6) Key for Printing Strains, Stresses and 

Strain Energies Per Unit Volume 
(if K£{{\) .NE. 0 - Print) 


If fiPLOT .EQ. 0 skip C^rd 6 and go to Card 7 


Card 6 {SEIO.S) 
Column 
PiO” PYSCL 
11-20 PZSCl. 
21-30 VNAX 
31-40 WMAX 

Plotter Control Card 
Contents 

Plot Y-Sca'ie Factor 
Plot Z-Scala Factor 
Maximum Allowable V-Displacement 
Maximum Allowable W-Displacement 

Card 7 {3(6X, D14.7)) 
Column 

Ej, Material Property Card 
^Contents 

1^.0 

27-40 

47-60 

C-bnstant term, A 
Linear tera, 8 
Quadratic tGrni» C 

Eu “ A • 

BT 5- CT^ 

Card 8 (3(6X, D14.7)) 

Ejfo ftotarlal Property Cerd 

Card 9 (3{6X, D14.7}) 

E-53 Material Property Card 

Ca.-d 10 (3{6X, D14»7)) 

G23 liaterial Property Ctird 

Ca»d 11 {3{6X, D14.7)) 

Gv> Material Property Card 

Card 12 (3(6X, 014.7)) 

Gj ^2 Material Proparty C^rd 

Card 13 (3(6X, D14.7)) 

NU23 Material Property Card 

Card 14 (3(6X, 014.7)) 

NUj 3 Material Property Card 

Card 15 (3(6X, D14.7)) 

NU;2 Material Property Card 

Card 16 (3(6X, 014.7)) 

Alpha 11 Material Property Card 

Card 17 (3{6X, D14.7)) 

Alpha 22 Material Property Card 

Card 18 (3{6X, 014.7)) 

Alpha 33 Material Property Card 

Card 19 (8010.5) 
Column 

•TT3" Ang(l) 

Angle Data Card 
Contents 

AnTfrenTo. i (In Degrees) 
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il-J5 Aog(2) 


Angle No, 2 (In Degrees) 


An g( NANG) Angle No. NANG (in Degrees) 


Card 20 (5X.15I5) 
Tclumn 


-rr 

6-10 

n-ib 

16-20 

21-25 

26-30 


B1 ank 

lANG(K) 

KOO(H.l) 

N00(N,2) 

N0D(N,3) 

N0D(N,4) 


Element Data Card(s) 

Cont ents 

{ Ca n firm ement Numbers in for Reference) 
Angle Number of Element N 
Nods #I of clement N 
ftode »Z of Element N 
Node #3 of Element N 
Node #4 of Element N 


(continue through NOD(N.NPE) for NPE > 4.) 


31-35 ISTRS(N,1) Stress and Strain Ojtput Code for Element N. 


0 

1 

2 

3 

4 

5 

C 

7 

8 
9 


None 
Side 1 
Side 2 
Side 3 
Side 4 
Center 
’•Node 1" 
"Node 2" 
"Node 3" 
"Nods 4" 


36-4G 

41-45 


ISTRS(N,2) 

ISTRS(N,3) 


71-75 ISTRS(N,9) 

♦Repeat Card 20 NEM Times 


Card 21 (5X/2D10.5) 
Column 

~yr5 Blank 
6-15 Y(N) 
16-25 Z(N) 


Nodal Data Card(s) 

Contents 

(Can l^jt S*ode Numbers in for Reference) 
Y-CoordinvTte of Node N 
Z-Coordinate ot Node N 


♦Repeat Card 21 NODS Times # 


♦If NSDF ,EQ. 0 - Skip Card 22 and Go on to Card 23 


Card 2 2 (215,010,5) 
licTlunin 
r~T~~ ND 
6-10 IDR 


SpecKied Nodal Degree cf Freedom Card(s) 
Contents^ 

Uo'3e Itiri’ber of Specifiec 0,0. F. 

Direction of Specified D.O.F. 

(I = u, 2 = Vj, and 3 “ w) 


11-20 VBOF(N) 


Spacifled Displacement Value 


^Repeat Card 22 NSDF Times 


* If KSSr .EQ. 0 Ski? Card 23 and Go on to Card 24 


i^rd 23 (215,010.5) 
('olumn 
"T^j f;0 
6-10 iDR 

n-20 VBSF(H) 


Specified Wodcl Force Card(s) 
Contents 

[idZfTl&'iFer of Specified Force 
Direction of ^ecified Force 
(1 “ u, 2 V, 3 * w) 

Specified Boundary Force Value 


^Repeat Card 23 NS5F Times 

*If NCON ,EQ. 0 - Skip Card 24 and Co on to Card 25 

Card 24 (215) Specified Constrained Displacensnt Card(c) 

'Tcriufr.n Contents 

KO T?oaT*1Iimbar of Specified Constrained 

Displacement 

6-iC lOR Direction cf Speti^ic'd Constrained 

Displaceiw^nt 

(I " u, 2 = V, and 3 « vj) 


* Repeat Card 24 iiCON Times 


Cdfd 25 (3(6X,D14.7)) 
Column 
'ITar' EPSX. 

27-40 TFREi 
47-60 TE!^ 


ISornal Strjin and Tetiperaturs Load Card 
Contents 

Ap‘pTTed'"Roi*mal, (X-'Jl recti on). Strain 
Stress Fre„ Temperatiira 
Analysis Temperature 



